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Gravity with extra dimensions and dark matter interpretation: 
Phenomenological example via Miyamoto-Nagai galaxy 

C. H. Coimbra-AraujcQ 

Instituto de Fisica Gleb Wataghin, Universidade Estadual de Campinas, Unicamp, 13083-970, Campinas, SP, Brazil. 

P. S. LetelieiQ 

Departamento de Matemdtica Aplicada, Instituto de Matemdtica, Estatistica e Computagao Cientifica, 
Universidade Estadual de Campinas, Unicamp, 13083-970, Campinas, SP, Brazil. 

A configuration whose density profile coincides with the Newtonian potential for spiral galaxies 
is constructed from a 4D isotropic metric plus extra dimensional components. A Miyamoto-Nagai 
ansatz is used to solve Einstein equations. The stable rotation curves of such system are computed 
and, without fitting techniques, we recover with accuracy the observational data for flat or not 
asymptotically flat galaxy rotation curves. The density profiles are reconstructed and compared to 
that obtained from the Newtonian potential. 
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I. INTRODUCTION 



q ■ There is strong observational evidence, primarily from dynamical and lensing effects, that galactic disks, cluster of 
q-i galaxies, and a smoothly distributed cosmological background point to the existence of the so called "dark matter". 

In galactic disks, where Newtonian gravitational theory would have been expected to be an excellent description, 
5jT)| accelerations of stars and gas, as estimated from Doppler velocities are much larger than those due to the Newtonian 
field generated by the visible matter in those systems (the plateau anomaly in rotation curves of galaxies) [H]. Rotation 
curves are the major tool for determining the distribution of mass in spiral galaxies, and are also important to study 
kinematics and to infer the evolutionary histories in galactic systems. Historically they are the most basic and classic 
■ manner to infer the presence of dark matter in galaxies (for a complete review about rotation curves, see for instance 
[2]). On the other hand, it is verified that cluster of galaxies are composed of three main components: ~ 5% in 
CN ■ mass is the optically luminous baryonic matter in hundred of bright galaxies; ~ 15% is in the form of a bright X-ray 
inter-cluster gas; and the remaining ~ 80% is some sort of non-baryonic "missing mass". Such techniques plus the 
. temperature fluctuations in the Cosmic Background Radiation have been regarded as confirming the dark matter 
' existence. 

a thin disk constructed from a space-time endowed with extra dimensions just provide the needed extra 
' parameters to construct, without dark matter, a configuration that mimics a generic and idealized axially symmetric 
J> \ galaxy. At this same reference, an impressive result was also obtained for gravitational lensing effects for a spherical 
cluster living within the same space-time. These outcomes are directly related to the dark matter problem, as it is 
explained above, in the sense that "dark matter" can be translated as a ignorance to explain why spiral galaxies have 
exotic rotation curves see e.g. 0, [f| and references therein, galaxy clusters have a greater amount of ray deflection 
than expected for gravitational lensing 0| (and also mass when confronted to virial theorem [l||), and the universe 
appears to have an unexpected few fraction of baryons produced during nucleosynthesis Q. 

In the present text it is given a phenomenological example where the previous thin disk model is naturally extended 
by using an isotropic configuration to construct a galaxy as mentioned in Q. Unlike the results obtained for the 
thin disk model, where the galaxy is interpreted as a flux and counterflux of geodesic particles, at the present article 
we obtain a stable configuration similar to a spiral galaxy, i.e., where we have a central bulge and a thick disk 
displayed in a particular density profile. In this way, the "artificial" galaxy constructed by the previous thin disk is 
substituted here by a richer model and a more realistic configuration. The design is totaly carried out by a AD isotropic 
configuration living in a multidimensional universe. It is constructed a 6D galaxy following the arguments about the 
simplicity of self-gravitating objects living in six dimensions (even number of space-time dimensions) presented in 
@|. The obtained Einstein equations allow the system to be solved by Miyamoto-Nagai solutions. Such solutions 
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are important because they describe a family of s elf-g ravitating configurations that can be seen as three-dimensional 
models for the distribution of mass in galaxies [1, flOl f 

The present work is organized as follows: in Section [ill the field equations are calculated from a AD isotropic metric 
plus extra terms and a Miyamoto-Nagai ansatz is used to solve them. In Section Hill a general equation for the circular 
geodesies in the planar part of the configuration is calculated and it is argued that such equation could approximately 
represent the rotation curves and in Section [IV] we use a general relativistic exact method to calculate the stability 
of such curves. Many spiral galaxies can be constructed from our achievements and in Section [V] it is presented an 
example for a well known galaxy (NGC 3198) and the density profile from the obtained rotation curves is recovered 
and compared with the Newtonian potential and other examples. Finally in Section IVIII we present some concluding 
remarks. In what follows we use c = 1 and G = 1 (and do not consider possible variations of c or G with space, time 
or number of dimensions) . 



II. FIELD EQUATIONS 



Consider a generalization where our universe has D = 4 + n dimensions. For an Einstein-Hilbert gravitational 
action we have 



s "lkJ dixdn y^- {4+n) 9 {A+n)R > (!) 

what leads to the field equations 

{i+n) G AB = -8tt {i+n) T AB , (2) 

where A,B = 0, 1, 4 + n — 1, y are the extradimensions and the indices (4 + n) tells about the multidimensional 
nature of the action. First of all, let us consider the case of axial-symmetric 4D space-times whose metric can be 
written in a isotropic form in cylindrical coordinates (t, i?, z, <p): 

ds 2 = e u(R,z) dt 2 _ e X(R, Z ) ((LR 2 + Az 2 + #2^2) (3) 

A general relativistic formulation for the Newtonian well known galaxy models can be written in the form of the 
Schwarzschild metric in isotropic coordinates (see e.g. [IB]). Also, as a form to extend the formalism developed in [jj, 
we introduce n extradimensional coordinates. For this case, it is showed in [E[ that non-exotic matter is only possible 
(including the fact about a satisfactory Huygens principle) if n is even. Let us consider e.g. the simplest case, where 
n = 2: 



(I-/) 2 
(1 + /) 2 

e" fc dx 2 - e fc dy 2 , (4) 



d,s 2 = ^^dt 2 -(l + /) 4 [di? 2 +dz 2 + i? 2 d^ 



where / = f(R, z) and k = k(R, z). The field equations (|A2|) yield the next expressions for the components of the 
energy-momentum tensor: 



with the constraint 



k 2 R - k\ = 0. (6) 

The general solutions for this last have the form 

k = k\{z — R) + k 2 or k = k\{z + R) + k 2 , (7) 

where fci and k 2 are constants, and for simplicity we will consider k 2 = 0. Also, with a such constraint, the AD part 
of pressure solutions are the same as obtained by Vogt and Letelier [Io| : 

T " = 4tt(1 + /)5(1 - /) { flzZ + 2f ' R ~ f ' z ' ' (8) 
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T ^ = T « = _ 4^(i^7)5(T 3 (//,«*- 3 /,«/,*) > ( 10 ) 

T ^ = 47r(l + /) 5 (l- f) ^ ^' RR + ^ ~ ^ R ~ ^ ' 



The extradimensional pressure part has the form 

Ml + /)s v j,rt " 1 J ' zz 1 i? 



= -t*. (is) 

The energy density is given by p = T t *, and the stresses (pressures or tensions) along a particular direction read 
Pi = —Tl when the energy-momentum tensor is diagonal. It is a surprising fact that the component T\ is proporcional 
to the usual Laplacian of the function / in flat cylindrical coordinates. Note that in the Newtonian limit when /< 1, 
Eq. |(5|) reduces to Poisson equation 

V 2 $=47T j O A r, (14) 

if the function / is related to the gravitational potential <& by 

$ 

/ = "2- ( 15 ) 

In this case, p —> pn and the energy conditions for the disk has no exotic matter (p + J^i^i > 0)- The energy- 
momentum tensor will be diagonal (T^ = Tjg = 0) provided / has the form 

C 

f= y/w(R)+g(zy (16) 

where C is a constant and w(R) and g(z) are arbitrary functions. Complementarily, and T£ will be equal (isotropic 
radial and azimuthal stresses) only if w(R) = R 2 . The density profile p can now be deduced, e.g., by Miyamoto-Nagai 
solutions 0, which represent stratifications of mass in the central bulges and in the disk parts of galaxies. In this 
case, the simplest gravitational potential that provides diagonal components is 

M 

*(R, *) = (17) 

y'i? 2 + (a + y/z 2 + fc 2 ) 2 

where a, b are positive constants. 

The corresponding three-dimensional density derived from Eq. (fl~4| is 



_ b 2 M aR 2 + (a + 3Vz 2 + 6 2 )(a + Vz 2 + fr 2 ) 2 
^ 4tt [i? 2 + (a + V2 2 + & 2 ) 2 ] 5 / 2 (z 2 + 6 2 ) 3 / 2 ) ' J 

and now the function f(R, z), according to Eq. lfl5|) . is 

M 

f(R, z) = (19) 
2tJr 2 + (a + Vz 2 + b 2 ) 2 

As we will see in next sections, the gravitational potential only can be approximated by f(R, z) in the Newtonian 
limit. By now, the true gravitational potential will be calculated from the circular velocity of some test particle in 
the system. 



4 



NGC 3198 (Begeman 1989) 
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PIG. 1: Rotation curves of NGC3198 implemented from gravitation with extra dimensions (GEDi, here with six dimensions) 
and a pure 47? Miyamoto-Nagai solution. The plots come only from the mass due to gas and stars, i.e., no dark matter 
is considered. The discrepancy at 5-7 kpc is due to the fact that our method is not a fitting technique and also the gas 
dynamics was not modeled. The result comes from circular geodesies of test particles around the designed configuration. We 
use b/a = 0.01 (that models a disky galaxy) and the stability parameter used here is fci = 10~ 6 . To address stability we use 
also (for extradimensional parameters) C x = 0.2 and C y = 0.8. 



III. APPROXIMATED ROTATION CURVES FROM CIRCULAR GEODESICS 

Now, the particles of such configuration describe trajectories that can be derived by calculating the system geodesic 
equations. Usually, when there are some stresses (and there are, both radial and azimuthal), the configuration can 
hardly be interpreted in terms of particles moving on circular geodesies. However, the assumption of geodesic motion 
is only valid for the case of a particle moving in a very diluted gas like the gas made of stars that models a galaxy 
disk. Given this important statement, it is possible to obtain the tangential velocities of the disk particles (i.e. the 
approximated planar rotation curves) from geodesic equations. Assuming R = and i = (the particles have no 
radial motion and for simplicity are confined in the surface z = 0), the metric (|4]) can be rewritten as 

jl^t 2 - (1 + /) 4 i?y - e- fc ± 2 - e k f = 1 (20) 

where x A = dx A /ds, which gives 

i 2 = (^-0 2 [i + (i + /) 4 i?V + c- fc i 2 + cV]. (2i) 

The Euler-Lagrange equations for x and y coordinates permit to calculate the following geodesic equations 

(e~ k xy = 0; e~ k x = C x , (22) 

(e k y)=0; e k y = C y , (23) 

where C x and C y are integration constants. We can fix the values of C x and C y by those calculated in [ll|, where it 
is obtained a stable planar configuration. 

Another equation is obtained by derivating l[2"Tj) in r and using eqs. (|22|) and lj23]) 



(1-/) , {I -If 



i 2 (24) 



(l + /)2 (l + /)3 
2R{1 + /) 3 [(1 + /) + 2Rf. B }^ 2 - k R (c- k x 2 - e k y 2 ) = 
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FIG. 2: Rotation curves of NGC3198. Also is considered a comparison between our model (a gravitation with extra dimensions, 
i.e. the GEDi Miyamoto-Nagai), a GEDi thin disk Q and other models or fittings [12, US, EH- The coordinates are the rotation 
velocity and the radius per mass. 



Eqs. f2"Tj) and ((24)1 form a system of equations whose variables are if 2 and t 2 . By solving the system, it is possible to 
calculate the rotation curves Vc by 



Va= I (25) 

9tt at y g tt t 2 



IV. STABILITY 



For such planar rotation curves, it is needed to calculate the stability. The most complete method is described in 
flTj |- In this case we compute the infinitesimal perturbations in the geodesic. The perturbation of the geodesic equation 
x A + Tg C x B x c = is done performing the transformation x A — > x A + A A — where A A = (St, Sr, Sip, Sz, Sx, Sy) are 
infinitesimal elements. As showed in it is possible to achieve to the following equations for perturbations: 

A A + 2T A c x B A c + T A CD A D x B x c = 0, (26) 

where Tg C are the Christoffel symbols and x A are proper time derivatives dx A /ds and can be written for a circular 
orbital motion as 

x A = (u t ,Q,Q,u t n,u t C x ,u t C y ), (27) 

where ft = Vc/R, Eq. <f25|l . Assuming oscillations in all directions with no vertical or extradimensional restrictions, 
we get 

A' 4 = (St, SR, Sz, Sip, Sx, Sy). (28) 
The non-null computed Christoffel symbols are: T\ R = T Rt , T\ z = r* t , rg, T% R , T* z = Tf R , I* , r*,, I* , 

pi? pz p2 p2 _ p2 p2 pz pz -pz p¥> _ p<£ T^ip T^<P 'px — T^ x p- T T^ x — py 

yy tt' RRi Rz 1 zRi 1 zzi W xxi yy> Rip ipR' zip 1 <pz i 1 Rx 1 xRi 1 zx ^ xzi Ry yR' 

T v zy = r^ z . Let x A be an equatorial circular geodesic in a stationary axisymmetric space-time {4]), i.e., the worldline 
x A = (t,R= const, ip = const + Qt,z = 0, x = const,?; = const). Substituting the four velocity l(27j) . and from the 
non-null computed Christoffel symbols, the components of Eq . (|26|1 and supposing that the solutions for St, SR, Sz, 
Sip, Sx and Sy have a form of harmonic oscillations, ~ e lKs , with a common proper angular frequency K , we find: 
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FIG. 3: (a) The Miyamoto-Nagai AD model Newtonian density pn in M Q pc -3 for NGC 3198 (the non-central disk, where the 
rotation curve anomaly appears), where R and z are given in kpc. Here k — K/u* . (b) The 3D plot of the recovered density 
profile for the present model. 



(St) + 2Y\ B u\5R) + 2T\ z u\5z) = 0, 

{5R) + 2T R u\St) + 
+[(r5, 1 



2r£ ¥ yn(M 



[(^tt,R "I" x 



^ 2 + r 



(<fc) + 2r£«*(<ft) 

+[(rti, a + + v xx z c 2 + r 

(<fy>) + 2r£ i? ft M t ( ( 5i?) + 2r£,nu*(<fe) = o, 
(Sx) + 2T x vR C x u\5R) + 2T x xz C x u t {8z) = 0, 
(6y) + 2Tl R C y u\5R) + 2Tl z C v u t {5z) = 0. 

The unique real solution for this homogenous system is 



t\2l 



R 

xx^R^x 



xx.R^x 



uu 



^(u^SR 



(29) 
(30) 

(31) 

(32) 
(33) 
(34) 



K 2 = ^ 



Rx 

3 



2 1 / 3 (-R 2 + R 2 ) 



2R\ - QRiR 2 + 27R 3 + y/A{-R 2 + 3i? 2 ) 3 + (2Rf - 9R^ 2 + 27R 3 ) 2 
1 



1/3 



(35) 



3(2V3) 



2R\ - 9i?ii? 2 + 27R 3 + \U{-R 2 + 3i? 2 ) 3 + [2R\ - 9R X R 2 + 27R 3 ) 2 



1/3 



where R\ 



{T R 



where T^b R u' 



AB,R 



r z AB > A u«, R 2 = (T R 



0(r 



AB,z l 



r R o 2 



r R r 2 

L xx,Ryx 



C 2 x C 2 {u l f, and 



u A u B ) and R 3 = 16T y yR T R y r- z T xx C 2 x C 2 ( 
Y R y i{Cy)(w*) 2 . The system is stable if the squared epicyclic 



AB,z 

(r^fl 

frequency k 2 = (K/u*) 2 is strictly positive. 

This occurs when the signal before the square root in Eq. Il35|) is positive and for k = k\(z + R) in Eq. (7]). More 
than this, the stability test may have an important role to determinate a stable range of values for the constants C x , 
C y and ki — this last a constant of integration which came from the presence of extradimensions for a non-exotic 
system, Eq. (7]). In the case of constants C x and C y , the values coincides exactly with the same obtained in [ll[. 
In the case of the constant ki, it must be tiny (ki ~ 10 -5 — 10 -7 ). The usual manner to determine parameters a 
and b (that are given in kpc) in Eq. lfl9|) is to obtain them phenomenologically by observing values of density for the 
galaxies, as done in [§]. Usually, the fraction b/a estimates how disky or spherical is the galaxy. Thus, parameters a 
and b can be rearranged in the new parameter b/a. Usually, disky galaxies have b/a ~ 0.01. Observation of densities 
plus the stability study present a complete manner to construct a galaxy for a system modeled from gravitation with 
extra dimensions. 



R kpc 



FIG. 4: The GEDi density contour plots for NGC 3198. We use b/a = 0.01 (that models a disky galaxy) and the stability 
parameter used here is fci = 10~ 6 . To address stability we use also C x = 0.2 and C y = 0.8. 



PROBING THE MODEL WITH REAL GALAXIES 



As showed in [8| one can model a real spiral galaxy by using a superposition of the central and outbound densities 
for different values of a and b. Here we will constrain as well the value of k as explained in the previous section. 
This is the same as write p = J2i P( a ^ I n other words, this is the same to write the function /, Eq. l[T9|) . by the 
superposition of J^i f( a i> For the case °f a spiral galaxy, say, it is suffice to write for the central bulge and for 
the disk galaxy part. It is possible indeed to calculate the rotation curves for many galaxies (including not only the 
spiral morphology). As an example, the spiral galaxy NGC 3198 can be modeled. Considering the surface density 
observed on the disk and the morphology characteristics [l2j and assuming that the bulge density is not so different 
than that observed for Milky Way H (where usually a is set to zero and b ~ lkpc; this means that b/a — > oo, or a 
total spherical ~ lkpc central bulge), we obtain a set of stable rotation curves (see Figs. [l]and[2|. 

The potential presented at Eq. (fl~7j) - i.e. a Miyamoto-Nagai ansatz - is purely Newtonian and was useful to 
calculate the form of the function / in the metric. Therefore it is necessary to calculate the true potential related to 
the problem. We can translate the potential form by calculating the density by 



V 2 $ 

where the potential is reconstructed from the planar circular geodesies (j25|) by 



$ = 



V 2 



(36) 



(37) 



where Vc is given in Eq. lf25]) . The density contour plots recovered from such equations is showed in Fig. [H and our 
object is very similar to a spiral galaxy. 

We can compared such result with other potentials known in literature. In e.g. it is listed some of possible 

galaxy potentials. A first one, used here as well (Miyamoto-Nagai potential-density pair) is compared to the density 
recovered from our circular velocities (see Fig. [5]). 

Another interesting potential to be compared is the Satoh's 



*%(R,z) 



M 



(38) 



ab 2 M 



Pm 



4ttS 3 (z 2 + b 2 ) 



R 2 + z* 
S 2 



(39) 



S=[R 2 + z 2 



+ a{a + 2^z 2 + b 2 )] 1 ' 2 . 



(40) 



Now, for same parameters a and 6, we plot those densities in Fig. [6l The bit difference between the reconstructed 
density profile and conventional Newtonian AD examples, shows that the extra dimensions affect remarkably the 
circular velocity profile but, although affect poorly some absolute values when compared to other potentials (what 
indeed does not disagree remarkably with observations), the density contours point to an object very similar to a 
spiral galaxy. 
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(a) (b) 
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FIG. 5: (a) The Miyamoto-Nagai AD model Newtonian density pn in Mqpc -3 for NGC 3198 (the non-central disk, where the 
rotation curve anomaly appears), where R and z are given in kpc. (b) The 3D plot of the recovered density profile for the 
present model. 
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FIG. 6: Comparison between the density recovered from our circular velocities (full line), a pure AD Myiamoto-Nagai profile 
(long dash line), and a Satoh profile (short dash line). Here the density is given in M© pc -3 and comes from the configuration 
slice where z = 0. 



VI. ROTATION CURVES ARE NOT FLAT! 

An essential observational remark is that real galaxy rotation curves are not flat. Thus the flat rotation curve 
paradigm should be dismissed for every galaxy model. Here it is possible to show that when the parameter ki is in 
the superior limit of stable curves, i.e., k\ ~ 10 -5 , the curves become non-planar. It is showed in Fig. [7l 

VII. CONCLUDING REMARKS 

Here we presented stable rotation curves of a configuration living in a 6D universe and the comparison with some 
potentials present in literature. Initially from a configuration whose density profile coincides with the Newtonian 
potential for spiral galaxies our model was constructed and a Miyamoto-Nagai ansatz was used to solve Einstein 
equations. The stable rotation curves of such system were computed and, without fitting techniques, we recovered 
with fidelity the observational data. 
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FIG. 7: The observed non-planar rotation curve of galaxy UGC 195 [T3|, the Miyamoto-Nagai 4D expected profile and the same 
when it is included more two extra dimensions in computation (the case of the present model). Parameters used: b/a = 0.01, 
C x = 0.2, C y = 0.9 and fci = 10" 5 . 

The present semi-phenomenological approach points to the evidence that a universe endowed with extra dimensions 
could explain the "missing mass" problem without a dark matter particle, at least for the case of rotation curves 
of galaxies, following the same results obtained in As explained in Q and also in [llj], although there is no 
fundamental theory presented, Eqs. lf2Tj) - (j23|) allow one to connect our model with some universal extra dimension 
like theory although here we relax about compactification. 

Extra dimensions arise here three main integration new constants (fci, C x and C v ), that are used as parameters of 
the model. The main goal is to find stable parameters using a geodesic perturbative approach. Thus we find that 
C x and C y should range as the conventional thin disk calculation obtained in [ll|, and ki is a very tiny parameter 
ranging from 1CT 7 — 10~ 5 . Miyamoto-Nagai configurations also introduce the parameter b/a and for disky galaxies 
we should use b/a ~ 0.01. 

Our model also, inside the range where curves are stable, reproduces the behavior of non flat galaxy rotation curves. 
The incredibly amount of theories that either imply or assume the existence of an observational scenario in which 
rotation curves of spirals are asymptotically flat, is clearly in contradiction with observational evidence. Thus, in this 
aspect, the present model explain planar and non-planar rotation curves. 

Here obviously there is no claim to replace the dark matter paradigm since the usual models of dark matter have been 
a relatively successful black box, fitting very well with many theoretical and observational issues: not only rotation 
curves and lensing, but also N-body numerical simulations, structure formation and the analysis of anisotropies of 
the CMB. To replace this black box by another one is not a simple matter, and requires its testing in a wide scope of 
theoretical and observational areas. 
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APPENDIX A: MOTION EQUATIONS FOR A TEST PARTICLE IN GRAVITY WITH EXTRA 

DIMENSIONS 

The Einstein-Hilbert gravitational action with more dimensions is given by 

S = [ d^d^V^+^g ( {A+n) R + C M ), (Al) 
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This leads to the field equations 

(i+n) G AB = {i+n) T AB , (A2) 

where A, B = 0, 1, 4 + n — 1, y are the extradimensions and the indices (4 + n) tells about the multidimensional 
nature of the action. From now, ' 4+ ™'Gab and ( 4+ ™'Tab will be simply called Gab and Tab (the same for the 
curvature tensor and scalar, and as well for the metric). 
The most general metric for the space-time given above is 

/ 9a0 | gab \ 

g(x a ) = , (A3) 

\ gap I g a b J 

where a, j3 = (0, ..,3) and a,b = (4, .., n), for any integer n ^ 4, and we are considering the conventional treatment to 
do the metric a function only of 3 + 1 coordinates. This metric, as written above, contains the explicit terms meaning 
the 3+1 universe and also the n terms plus crossed components. In fact, Eq. (| A3[) can be as well rewritten, for 
convenience, in a different way as 

gAB = g a p$ A V B + 9abS A 8 B + gab^A^B + gapS A S B , (A4) 

for A, B = (0, .., 3 + n) and where dj are the conventional Kronecker symbols. 

Considering the particular case where the metric is diagonal, it is possible to find the derivatives and consequently 
the curvature terms (see [i| for the explicit calculation) . The equations of motion for such system are calculated as 



+ {<&} = ^gab n g^N c g ac N d g M , (AS) 

where N a means a vector that contains, e.g., the parameters C x and C y (for the case of 6D). The metric elements 
should be calculated by a new Poisson equation (plus boundary conditions and initial values) that arises from the 
new terms in Einstein equations. The Myiamoto-Nagai example in the present paper is a particular case for such 
approach. All the discussion about those terms are presented at [§]. 
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